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Abstract 

In our previous paper [7j we have shown that the theory of circulant ma- 
trices allows to recover the result that there exists p + 1 Mutually Unbiased 
Bases in dimension p, p being an arbitrary prime number. Two orthonormal 
bases £>, B' of C d are said mutually unbiased if V6 G B, W G B' one has that 

\b-b'\ = -±= 

Vd 

(b ■ b' hermitian scalar product in C d ). In this paper we show that the theory 
of block-circulant matrices with circulant blocks allows to show very simply 
the known result that if d = p n (pa prime number, n any integer) there exists 
d + 1 mutually Unbiased Bases in C d . Our result relies heavily on an idea 
of Klimov, Muhoz, Romero As a subproduct we recover properties of 

quadratic Weil sums for p > 3, which generalizes the fact that in the prime 
case the quadratic Gauss sums properties follow from our results. 



1 INTRODUCTION 

The Mutually Unbiased Bases in dimension d are a set {£>i, B^} of orthonormal 
bases in C d such that for any G b^ ' G £>&/ one has 



bf ■ 



-j=, Vj,f = l,...,d } W^k = l,...,N 



where b ■ b' = Ylj=i ^ s ^ ne USU& 1 scalar product in C d . 

This notion of mutually unbiased bases emerged in the seminal work of Schwinger 
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[14| and turned out to be a cornerstone in the theory of quantum information. Fur- 
themore it is strongly linked with the theory of Hadamard matrices |9j and to the 
Gauss Sums properties. 

An important problem is the maximum number of mutually unbiased bases 
(MUB) in dimension d. The problem has been completely solved for d = p n where p 
is a prime number, and n any integer, in which case one can find N = d+1 MUB's 

In a previous paper [7] (hereafter refered to as I) we have provided a construc- 
tion of d + 1 MUB's for d a prime number using a new method involving circulant 
matrices. Then the MUB problem reduces to exhibit a circulant matrix C which is a 
unitary Hadamard matrix, such that its powers are also circulant unitary Hadamard 
matrices. Then using the Discrete Fourier Transform which diagonalizes all cir- 
culant matrices, we have shown that a MUB in that case is just provided by the 
set of column vectors of the set of matrices {F^, 11, C, C 2 , C d_1 }. Properties of 
quadratic Gauss sums follow as a by-product of the method. 

The present paper is a continuation of I, where we consider d = p n . Here 
circulant matrices are replaced by a set of block-circulant with circulant blocks 
matrices. Again the discrete Fourier transform which in this case will be simply 
F = F p eg) F p ® ... (g) F p will play a central role since it diagonalizes all block-circulant 
with circulant blocks matrices. We follow an idea of [UJ to define them. The new 
result developed here is that these block-circulant matrices with circulant blocks 
together with F will solve the MUB problem in that case. 

Let Bk = j&o^i ^ e orthonormal bases. Then in any given base, 

they are represented by unitary matrices B k . Taking B to be the natural base, we 
have that 

b 3 k ■ by } = (BlBk')jj' 

Thus in order that the bases Bk be unbiased, we just need that all the unitary 
matrices B^Bk', k ^ k' have matrix elements of modulus dr x l 2 . Such matrices are 
known as unitary Hadamard Matrices ([9]). 

2 THE SQUARE OF A PRIME 

Let p be a prime number. One defines a primitive p-th root of unity : 
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The Discrete Fourier Transform in C p is 
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Definition 2.1 Consider a d-periodic sequence a , a±, aj-i, a , a±, .... 

(i) A d x d matrix D is diagonal and called diag(o , ...,0,4-1) if its matrix elements 

satisfy 

D j:k = a k 5 jtk , Vj, k = 0, 1, d - 1 

^ A d x d matrix C is called circulant and denoted C = circ(ao, aa-i) if its 
matrix elements satisfy 

Cj : k = d{d~l)j+k 

Thus it can be written as : 

I a a 1 a 2 ■ ■ a d A 

q _ a-d-i clo a>i ■ ■ a>d-2 

\ ai a 2 a 3 . . a / 

(mj ^4 diagonal and circulant matrix must be a multiple of the identity matrix 1. 
(wj A d 2 x d 2 matrix is said to be block- circulant if it is of the form 

C = circ(Co,Ci,...,C d _i) 

where the Cj are d x d matrices. 

(v) It is block- circulant with circulant blocks if furthermore the Cj are circulant. 

(vi) A d x d matrix H is a unitary Hadamard matrix if 

d-i 

\H hk \=d- 1 ' 2 , and Y, H ik H k,i = Sj,i 

We define the following p x p unitary matrices 

X = circ(0,0,...,l) 

Z = diag^w,...,^ -1 ) 
They obey the u>- commutation rule : 

Lemma 2.2 (i) X p = Z p = 1 

(ii) 

ZX = ujXZ 

(Hi) Furthemore one has 

f p xf; = z 



(i) and (ii) are obvious. For a proof of (iii) see [Sj. 

Proposition 2.3 Let C = circ(ao, 

(i) One has 

p-i 

C = Y,a k XV- k 

k=0 

(ii) The discrete Fourier transform diagonalizes the circulant matrices : 

F p CF* = diag(a , ...,a p _i) 

where 

p-i 

k=0 

(iii) The set of circulant p x p matrices is a commutative algebra. 
Proof : 

p-i 

f p cf; = J2*kZ p ~ k 

k=0 

But Z~ k = diag(l, ui~ k , ...,cu~ k ( p ~ 1 ^), hence (ii) follows, (iii) is a consequence of (i). 

Corollary 2.4 If the sequence {a>k}ke¥ d ^ s suc ^ ^hat 

\dk\ = d~ l t 2 and \dk\ = 1 

then C is a circulant unitary Hadamard matrix. 

Proof : C is unitarily equivalent to an unitary matrix if |%| = 1. Furthemore 
\ a k\ = dr x l 2 , hence the result follows (see j7]). 

2 

The discrete Fourier transform in C p is defined as follows 

F = F P ®F P (2.1) 

It has the following important property (similar to the property that the discrete 
Fourier transform diagonalizes all circulant matrices) : 

Proposition 2.5 (i) F is an unitary Hadamard matrix. 

(ii) All block- circulant matrices C with circulant blocks are diagonalized by F : 

FCF* = D 
where D is a p 2 x p 2 diagonal matrix. 
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For a proof of this result see [H]- 

We shall be interested in finding block-circulant with circulant blocks unitary 
matrices in C p that are Hadamard matrices. An example is of course C ®C where 
C, C are unitary circulant Hadamard matrices. 

For p a prime number, denote by F p the field of residues modulo p. The correspond- 
ing Galois field GF(p 2 ) is defined as follows. For any p there exists an irreducible 
polynomial of degree two, with coefficients in ¥ p so that if we denote by a a root of 
this polynomial, 

GF(p 2 ) = {ma + n} m ^ ¥p 

The product 9 ■ 9' G GF{p 2 ) for 9, 9' G GF{p 2 ) is obtained using the irreducible 

polynomial which expresses a 2 in terms of a and 1. 

The additive characters i n GF(p 2 ) are defined as follows: 

Definition 2.6 The additive characters on GF(p 2 ) are : 

/n\ ( 227T 

X{9) = exp ( — — trfe* 

where 



p 

tr9 = 9 + 9 P 



Lemma 2.7 They satisfy 

(i) x (e + 9>) = x (0)x(0') 

(ii) One has : 



(Hi) 



8eGF(p 2 ) 

£ x(o-e')=P%,o (2.2) 

0'eGF( P 2 ) 

We take as natural basis in C p2 the set of states labelled by 9 G GF(p 2 ), in the 
following order : 

B = {|0),|a),|2a),...,|l),|l + a),...,|p-l),...,|p-l + (p-l)a)} 
Labelled by 9 G GF{p 2 ) we define a set of unitary operators in C p2 such that : 
Definition 2.8 (i) The set of operators T = {Xg} eeGF ( p2 j obeys 

X e \9') = \9 + 9'), W G GF(p 2 ) (2.3) 
(ii) The set of diagonal operators T' = {Ze} e&GF ^ p2 -j obeys 

Z \9') = x{e-9')\9) ) WeGF(p 2 ) (2.4) 
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They obey : 
Proposition 2.9 (i) 

Z e X e , = x(0-0')X e ,Z e (2.5) 
(ii) The operators in J 7 , T' obey the group commutative property : 

Xo+e> = XqXqi = Xg/Xg, Zq + qi = ZgZgi = Zg/Zg (2-6) 

(Hi) Z = X = 1 

Proof : Take any ip G GF(p 2 ). Then 

Z X ,\<p) = X (0 ■ (0' + y))\9' + ip) = x (9 ■ 6')X e ,Z e \ v ) = x (9 ■ 9') X (9 ■ <p)\9' + <p) 
We also have (ii) : 

ZeMv) = X({0 + 0') ■ <p)\<p) = X(0 ■ <p)x(& ■ <f)W) = ZoZo>\<p) 

□ 

Theorem 2.10 (i) 

?={x m ® xn m , neFp 

More precisely one has 

X ma+n = X n ®X m (2.7) 

(ii) All operators in T are represented by unitary block- circulant with circulant blocks 

p 2 x p 2 matrices. 

(Hi) 

r = {z m ® z*} mineWp 

(iv) All operators in T 1 are represented by diagonal p 2 x p 2 matrices. 

(v) For any 9' e GF(p 2 ) there exists a 9 £ GF(p 2 ) such that 

FX e >F* = Z e (2.8) 

Proof of (i) : It is enough to see that X a — 1 ® X and X\ — X eg) 1 since the 
other matrices X ma+n will be given by the chain rule : 

ma+n — -™- a 

But these are obviously block-circulant with circulant blocks matrices. 
One has for 9' = ma + n : 

FX ma+n F* = FX™X^F* = (F p ® F p )(X n ® X m )(F* p <g> F* p ) = Z n <g> Z m 

which is Zq for some 9 e GF(p 2 ). □ 
One recalls a famous result [8] : 
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Proposition 2.11 All block- circulant matrices with circulant blocks commute and 
are diagonalized by F . 

Proof : It follows from (12.81) that if C = Y^,g> Xg>Xgi is a block-circulant matrix with 
circulant blocks, one has 

FCF* = XqiFXqiF* = ^ XgrZf^gij 

e>&GF(p 2 ) e>eGF{ P 2 ) 

which is a diagonal matrix. 

To find the MUB's in dimension p 2 it is enough to exhibit a partition of the set 
of unitary operators : 

£ = {ZgX e/ } efi/eGF{p2) 

into a set of commutant families : We define 

= A W 

One wants : 

£ = Jo |J Cg U {11} 

9£GF(p 2 ) 

The family Cg will be defined as follows : 
Definition 2.12 Let for any 9 G GF(p 2 ) 

£g = {Zg,Xg.g,} e , eGF{p2) 

Define 

Cg = £g\{l} 

Proposition 2.13 (i) £ = T' 

(ii) Eg is a commuting family W9 G GF{p 2 ). 

(Hi) £ = Tq \Jg E GF(p 2 ) U {1} is a partition of 8. 

Proof : (i) is obvious. 

(ii) W, 9" G GF(p 2 ) one has 

Zg> Xg.gr Zgr'Xgn.g = x{~^ ' $ ' $ )^0'+0"^0-(0'+0") = ZgnXg.gn ZgiXg.gn 

(iii) and JF contain p — 1 elements. The classes Cg for different 9 are disjoint. 
Therefore 

U & 

9&GF{p 2 ) 

contains p(p — 1) elements. One has : 

p — 1 + p(p — 1) + 1 = p 2 
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which is the total number of elements in 8. 

Since all the unitary operators in Cg commute, they can be diagonalized by the same 
operator Rg. In the above cited work [11] they are defined as "rotation operators". 
In fact we shall see that they are represented in the basis B by block-circulant with 
circulant block matrices. The first main result of this paper is the following : 

Theorem 2.14 (i) There exists a set {Re} d( z GF ^ of unitary operators which di- 
agonalize all the operators of the class Cg, V# E GF(p 2 ). 

(ii) The operators Rg for 9 ^ are represented in the basis B by block-circulant with 
circulant block matrices which are unitary Hadamard matrices. 
(Hi) For p > 3 they obey the group law : 

R e+e/ = R e R e ,, V0, 9' E GF(p 2 ) 

Proof : It is enough to show that for any 9 E GF(p 2 )\ {0} the Rg can be expanded 

as 

Re= X e> X °' ( 2 - 9 ) 

0'eGF( P 2 ) 

since they will automatically represented in the basis B by block-circulant with 
circulant blocks matrices. One has to check that 

Rg 1 Zg/Xg.giRg = /Jg^g/Zg/, \/9' E GF(p 2 ) 

Since the operators ZgiXg.gi are unitary, the [ig^ are necessarily complex numbers 
of modulus one. But 

Zg'Xg.gi XffiXg" = Zgi Xg), ' Xg" + g.g> = [Iggt Xgl Xgm Zgi = flg^'Zgi x(~ 9' •9"')Xgn, Xgi, 

gn gu g,„ gm 

Equating the coefficients of Xgm in both sides we get 

^0'»-0.0' = ^e,e'X{—Q' ■ 9"')X g e ,l 
Taking 9"' = and assuming that Xq = p^ 1 , V9 E GF(p 2 ) we get 

A -g.g> — P ^6,0' 

or equivalently, since 9^0 

A gi — P fJ-e-e- 1 ^ 1 
This proves that all the A^ must be of modulus p~ x . 

Therefore since all the Xg are represented by unitary matrices that have non-zero 
elements (actually 1) where all the others have zeros, and since every matrix element 
of Rg is of the form A^ for some 9' E GF(p 2 ) } this proves that all the Rg are 
represented by Hadamard matrices. 
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Now we have to check the compatibility condition. We reexpress it in terms of ^e,e'- 
Supressing the index 9 in the fj,g t g/ for simplicity, we need to have W, 9" G GF(p 2 ) 

fJ--e-i(0"-96') = ^6'H-e- 1 e"X(—@' ' 9 ) 

or in other terms 

fie'+e" = H$'fie>'X(Q ■ & ■ 9") (2.10) 
But this results easily from the group property of the Xg's and Z 9 's ( 12.31 12.4[) : 

Rg 1 Zgi + e» Xo.(Qi + 0ii)Rg = x(99'9")R e 1 Zg/XggiRgR d 1 ZgnXggnRg = x(99' 9") flg^' flg t g/i Zg, 

In [TT] it is shown that for p > 3 the solution of ( 12.101) with /i<9 = 1 is 

H6,v = Xp~ 1 0-& 2 ) (2-11) 

Thus we deduce that 

\ e ) = p - 1 X {2- 1 9- 1 -{9'f) (2.12) 

We now prove the unitarity of Rg. For 9 = this is obvious since Rq = t. It is 
enough to check that for 9^0 one has : 

2 , \ x w ) A e'+e" — °e",o 

9>eGF(p 2 ) 

One has : 



/ j y A 8' ) A e>+e>> — 2 

6»'eGF(p 2 ) ^ 0'eGF(p 2 ) 



l 1 6-O- 1 -O'l J '9-0- 1 -(O'+O") 



= — ^2 \f 1 d-9- 1 -e'\ 2 ^e,-e- 1 -e"X(9 1 9'9") = \i e _g-x,g„ — ^ = <W,o 

P 6»'eGF(p 2 ) P 9>£GF(p 2 ) 

where we have used (I2.10p and (12. 2p . together with the fact that 9~ 1 ■ 9" = implies 
9" = due to the field property of GF(p 2 ). Thus one has 

r- 1 = r*, v0 e GF( P 2 ) 



(iii) The group law Rg + g/ = RgRg> for p > 3 has been established in [TT]. For p = 2 
it has to be suitably modified as shown in [TT] (see remark below). Let us see how 
it works for p > 3 : 

(Rl) n (R* a ) m Z g ,X 

(n+ma)- 

giR^R™ — (R\) n (R* a ) m ZgiXgi. raa R!^X n giR r l — fi(ma,9')(Rl) n Zg'X ni 

— ^ma,9'^n,9'Zgi = fJr ma +n,9' % 9' 

holds for p > 3 since 

Hma+n,9' = ^ma,9' Hn,9' , V0' G GF(p 2 ) 

which easily follows from (12.111) for p > 3, and the additivity of the characters. 
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Remark 2.15 For p = 2 the group law is not satisfied (see [11])- One has instead 
a very similar property (modified group law ) : 

R a R a+ i = R a+ iR a = R\ 

R a R\ — Ra+lX a+ \ 
Ra+\R\ — RaX a 

R 2 a — x a+ i 

p2 _ y 
R\ = X\ 

The second main result of this paper is the following : 

Theorem 2.16 The set of operators {F, Re}g e cF(p 2 ) defines a set of p 2 + 1 MUB's 
in C p \ 

Proof : Each Rg is represented by an unitary Hadamard matrix, and so is F. 
Due to the group property, R* e Rg> = Re>-e thus is an unitary Hadamard matrix. It 
remains to show that RgF* is an unitary Hadamard matrix V# £ GF(p 2 ). But we 
have 

R e F* = F*D 9 , y9 E GF(p 2 ) 

with Dg an unitary diagonal matrix, since F diagonalizes all block-circulant matrices 
with circulant blocks. The product of Dg with the unitary Hadamard matrix F* is 
obviously an unitary Hadamard matrix. 

In the case of dimension 2 2 one has instead of the group property that 

^9, e', 36" such that R* e R , = Rg^gXg,, 

Since Xgn has exactly one non-vanishing element ( 1) on each line and column, the 
product Rgi-gXgn is indeed an unitary Hadamard matrix. □ 



3 THE CASE d = p n 

The case d = p n with general power n can be treated similarly. There is a notion of 
block-circulant matrices with block-circulant blocks and building block pxp matrices 
which are circulant which generalizes the case n = 2. These are diagonalized by the 
Discrete Fourier Transform F in C p ™ which is 

F = F p ® F p ® ... <g> F p 

(n times) which is obviously an unitary Hadamard matrix. 

The Galois field GF(p n ) is defined through the irreducible polynomial which is of 
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power n, with coefficients in F p . Thus elements of the Galois field GF(p n ) are of 
the form 

n-1 

e = Cid 

where q G F p and a is a root of the characteristic polynomial. 
The characters are 

X{0) = exp (jj-tr(fi)\ 

where 

tr(G) =6 + 6 p + ... + 9 pn ~ 1 

For any 9 G GF(p n ) the operators R 9 are defined in [11] and the Rg obey 

a group law (resp. a modified group law) if p > 3 (resp. p = 2). 

As previously we have X e G {X fcl <g> X fc2 ® ... ® X kn , with % G F p } and 

z e e{z^®Z k *®...®z^} ki&p 

and 

8'&GF(p n ) 

with Ag, of modulus The Ag/ are given for p > 3 similarly to (12.121) by 

\ ( e 6 ) =p- n '\{2- 1 e-\e') 2 ) (3.13) 

All the results of the previous section are easily generalized. 

Theorem 3.1 The unitary Hadamard matrices F,Rg for 9 G GF(p n ) define a set 
ofp n + l MUB's m C pn . 



4 Weil sums for d = p n , p > 3 

The Weil sums in dimension p n are the equivalent of the Gauss sums for d = p (p 
prime number). The characters x($)> £ GF(p n ) replace the powers uj n , n G F p . 
Usually in the literature (see for example [H]), the Weil sums properties are used 
to solve the MUB problem. Here, as in [7], we do the converse. In the previous 
sections we have constructed the d+ 1 bases, and we shall deduce the Weil sums 
properties from this construction. 



Theorem 4.1 Let p > 3. Then for any 6 G GF(p n ) \ {0} and any 9' G GF(p n ) we 
have 



6"eGF(p n ) 



(4.14) 
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Proof : The matrix elements of a given row in F are of the form 

1 ?X(0' ■ 0") 



for some 9' E GF{p n ). Since 

Re = ^= E Xim- 1 ■ (9") 2 )X e „ 
e"eGF( P n ) 

the matrix elements of the first column of the matrices Re, 9^0 are of the form 

1 



X{2 9 1 • (9") 2 )e"eGF( P ") 



where we have used (13.131) . But the elements (20) _1 when 9 e GF(p n )\{0} span 
all 9\ G GF(p n )\ {0}. We have established that the matrix FRg is Hadamard. This 
implies that all its matrix elements have modulus p~ n l 2 . All matrix elements of the 
first column of FRq are thus of the form 

P- H E X(9i ■ (9"f + 9' ■ 9") 

6»"eGF(p") 

with 9\ = (29)~ 1 . Writing that its modulus is p~ n / 2 yields equ. (I4.14p . 
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